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ABSTRACT 

We present the results of our broadband spectral analysis of 42 SGR J1550— 5418 bursts simultaneously detected 
with the Swift /X- ray Telescope (XRT) and the Fermi /Gam ma-ray Burst Monitor (GBM), during the 2009 January 
active episode of the source. The unique spectral and temporal capabilities of the XRT windowed timing mode have 
allowed us to extend the GBM spectral coverage for these events down to the X-ray domain (0.5-10 keV). Our 
earlier analysis of the GBM data found that the SGR J1550— 5418 burst spectra were described equally well with 
either a Comptonized model or with two blackbody functions; the two models were statistically indistinguishable. 
Our new broadband (0.5-200 keV) spectral fits show that, on average, the burst spectra are better described with two 
blackbody functions than with the Comptonized model. Thus, our joint XRT-GBM analysis clearly shows for the 
first time that the SGR J1550— 5418 burst spectra might naturally be expected to exhibit a more truly thermalized 
character, such as a two-blackbody or even a multi-blackbody signal. Using the Swift and RXTE timing ephemeris 
for SGR J 1550— 5418 we construct the distribution of the XRT burst counts with spin phase and find that it is not 
correlated with the persistent X-ray emission pulse phase from SGR J 1550— 5418. These results indicate that the 
burst emitting sites on the neutron star need not to be co-located with hot spots emitting the bulk of the persistent 
X-ray emission. Finally, we show that there is a significant pulse phase dependence of the XRT burst counts, likely 
demonstrating that the surface magnetic field of SGR J1550— 5418 is not uniform over the emission zones, since it 
is anticipated that regions with stronger surface magnetic field could trigger bursts more efficiently. 

Key words: pulsars: individual (SGR J1550— 5418, IE 1547.0-5408, PSR J1550— 5418) - stars: neutron - 
X-rays: bursts 

Online-only material: color figure 


1. INTRODUCTION 

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs) are the observational manifestations of magnetars — 
isolated neutron stars possessing extreme magnetic fields, B 
> 10 14 G (Duncan & Thompson 1992; Kouveliotou et al. 
1998). Besides being bright X-ray sources, SGRs and AXPs 
emit intense bursts in hard X-rays/soft y-rays on a highly 
unpredictable frequency with peak luminosities ranging from 
10 38 erg s 1 to >10 47 erg s 1 . These energetic events are 
attributed to the cracking of the solid neutron star crust by 
magnetic stress buildup (Thompson & Duncan 1995) or to 
magnetic field line reconnection (Lyutikov 2003). For detailed 
reviews on SGRs and AXPs, see Woods & Thompson (2006) 
and Mereghetti (2008). 

SGR J1550— 5418 was discovered as a point source with the 
Einstein Observatory while searching for X-ray emission from 
radio emitting supernova remnants (Lamb & Markert 1981). 
The source was suggested to be a magnetar candidate by the 


1 1 Current address; The Open University of Israel, 1 University Road, POB 
808, Ra'anana 43537, Israel. 


similarity of its persistent X-ray spectrum to AXPs and its asso- 
ciation with a young supernova remnant (Gelfand & Gaensler 
2007). Its magnetar nature was confirmed with the detection of 
radio pulsations with P = 2.096 s and P = 2.318 x 10 -11 , 
corresponding to an inferred surface dipole magnetic field 
strength of 2.2 x 10 14 G (AXP 1E1547.0 - 5408; Camilo et al. 
2007). An accurate source location was also derived from the 
radio image, (J2000) R.A. = 15 h 50 m 54)ll ± (T0L deck = 
—54° 18'23"7 ± O'.' 1 (Camilo et al. 2007). X-ray pulsations at the 
same spin period were later found with a deeper XMM-Newton 
observation (Halpern et al. 2008). No bursts were detected 
from SGR J1550— 5418 until 2008 October, when both the 
Swi/i/Burst Alert Telescope (BAT) and the ft'ivni/Gamma- 
ray Burst Monitor (GBM) were triggered by numerous bursts 
from the source (Israel et al. 2010; von Kienlin et al. 2012). 
SGR J1550— 5418 entered an episode of more active burst- 
ing in late 2009 January, and no more bursts were detected 
after 2009 April. During these active episodes, several high 
energy instruments, such as the Swift/ BAT and X-ray Tele- 
scope (XRT), the Fermi/ GBM, the Rossi X-ray Timing Explorer 
(RXTE)/Proportional Counter Array, and the International 
Gamma-Ray Astrophysics Laboratory (INTEGRAL) /Imager on 


1 


The Astrophysical Journal, 756:54 (12pp), 2012 September 1 


Lin et al. 


board the INTEGRAL Satellite (IBIS) and the SPectrometer 
on INTEGRAL (SPI) recorded hundreds of bursts (Mereghetti 
et al. 2009; Savchenko et al. 2010; Kaneko et al. 2010; Scholz 
& Kaspi 2011; van der Horst et al. 2012; von Kienlin et al. 
2012). Following these burst active periods, both the persistent 
X-ray emission characteristics and the spin-down behavior of 
the source changed remarkably (Enoto et al. 2010; Ng et al. 
2011; Bernardini et al. 2011; Scholz & Kaspi 2011; Dib et al. 
2012; Kuiper et al. 2012). 

The spectral properties of SGR J1550— 5418 bursts have 
been extensively studied using individual instruments: SPI 
(Mereghetti et al. 2009), IBIS (Savchenko et al. 2010), BAT 
(Israel et al. 2010), XRT (Scholz & Kaspi 2011), and GBM 
(van der Horst et al. 2012). The XRT data cover a relatively 
narrow energy range (0.5-10 keV) of the spectrum of a typical 
SGR burst. Nevertheless, Scholz & Kaspi (2011) modeled 
SGR J1550— 5418 burst spectra using the XRT data only in the 
energy range of 0.5-10 keV with a single power law and found 
an average photon index of 0.17 ± 0.33. They also reported 
that there is a slight anti-correlation between the photon index 
and the absorbed X-ray flux. BAT provides a spectral energy 
coverage for SGR bursts from 15 keV to 150 keV. Israel et al. 
(2010) found that the spectra of BAT detected bursts can be well 
described by a single blackbody function with temperatures 
~10 keV. Finally, in van der Horst et al. (2012) we derived 
the spectra for a large set of SGR J1550— 5418 bursts detected 
with GBM in 2009 January using several continuum models. 
We found that in a slightly broader energy range (8-200 keV), a 
Comptonization model or the sum of two blackbody functions 
(BB+BB) can fit the SGR J1550— 5418 burst spectra equally 
well. Note that these two models were also used to describe the 
spectra of other magnetar bursts in a similar energy range and 
revealed intriguing physical insights into the burst phenomena 
(Feroci et al. 2004; Olive et al. 2004; Israel et al. 2008; Lin et al. 
2011). 

In this paper, we combine the spectral data of bursts observed 
simultaneously with XRT and GBM to investigate their spectral 
characteristics over a broader energy band (0.5-200 keV). In 
particular, we concentrate on the two most plausible representa- 
tions of SGR burst spectra, namely the Comptonization model 
and the BB+BB. Focusing on bursts with data collected over 
broader spectral bands enhances the chance to discriminate be- 
tween different spectral models. In Section 2, we describe both 
the XRT and GBM observations of SGR J1550— 5418 bursts 
and the selection of their common events’ sample. We present 
the data reduction and analysis in Section 2.3. The broadband 
spectral analysis results and their physical interpretation are 
presented in Sections 3 and 4, respectively. 

2. OBSERVATIONS AND SIMULTANEOUS 
BURST IDENTIFICATION 

2.1. Observations 

Swift/XRT is an X-ray imaging spectrometer sensitive to 
photons in the 0.3-10 keV energy range (Burrows et al. 2005). 
The telescope is operated either in photon counting (PC) mode or 
in windowed timing (WT) mode. Both modes provide the same 
spectral capabilities, however, the temporal resolution of the PC 
mode data is about 2.5 s, too coarse for the very short (<100 ms) 
SGR bursts. Therefore, we only employ here XRT observations 
performed in WT mode because of its 1.7 ms readout time 
well suited to the SGR burst durations. XRT monitored the 
source with 46 pointed observations in WT mode between 2008 


Table 1 

XRT Observations of SGR J1550— 5418 with Simultaneous Events with GBM 


Observation ID 

Date 

Start Time 

Exposure 2 

(ks) 

00340573000 

2009 Jan 22 

02:26:22 

6.38 

00340573001 

2009 Jan 22 

09:18:28 

9.45 

00030956035 

2009 Jan 30 

17:49:33 

2.97 


Note. a In WT mode. 


October 1 and 2009 April 30, covering burst active episodes 
of SGR J1550— 5418. These XRT observations were densely 
concentrated around the most burst active period (2009 January) 
with a total exposure time of ~175 ks. Three out of these 
46 pointings (listed in Table 1) included bursts simultaneously 
detected with GBM, as described in detail in the next section. 

The Fermi/ GBM monitors the entire sky (excluding the 
portion occulted by the Earth) in the energy range from 8 keV 
to 1 MeV with 12 Nai detectors and in the 0.2-40 MeV energy 
band with two BGO detectors. In its trigger mode, GBM records 
time-tagged event (TTE) data with high temporal and spectral 
resolution of 2 /is and 128 energy channels, respectively. The 
trigger readout lasts for 600 s (see Meegan et al. 2009 for 
more details of the instrument and data types). Using the 
same burst finding algorithm described in van der Horst et al. 
(2012), we searched for triggered and untriggered events from 
SGR J 1550— 5418 during the active periods of 2008 and 2009. 12 
In total, we identified 692 bursts out of which, 458 events had 
TTE data. We only used Na i detector TTE data for our spectral 
analyses, as typical SGR bursts are not detected above 200 keV, 
and last only for a fraction of a second. 

2.2. Identification of Simultaneous Events 

To identify the events observed simultaneously with GBM 
and XRT, we compared the times of the 458 GBM bursts with 
the time intervals of the XRT observations in WT mode and 
found 87 common bursts. Note that some SGR bursts have 
multi-peaked time profiles and each peak was labeled as a burst 
in our initial untriggered event search. We used the convention 
described in van der Horst et al. (2012) to determine whether 
multiple peaks constituted a single event — namely, we requested 
that the time difference between successive burst peaks was 
less than a quarter of the spin period of SGR J1550— 5418 
(~0.5 s). As a result, we obtained 66 SGR J1550— 5418 events 
simultaneously observed with GBM and XRT. For each burst, 
we plotted the burst light curves seen with GBM and XRT to 
determine the time interval that includes the main emission 
episode and is used as our spectral extraction interval (see 
Figure 1). We excluded nine dim bursts that had less than 
40 counts in the XRT data during the burst interval, statistically 
not enough for spectral analysis. Additionally, we excluded five 
very bright bursts that saturated the high-speed science data bus 
of GBM (see also van der Horst et al. 2012). 

We checked all common bursts in the XRT data for pile- 
up. We regenerated the level 2 data by including photons 
in all grades (0-15) from level 1 data at the position of 
SGR J1550— 5418, using the standard XRT data processing 
tasks in HEASOFT. Then we calculated the average readout time 
for each detection area. 13 We found that, besides the five GBM 


12 The list of GBM triggered events is available from http://gammaray.nsstc. 
nasa.gov/gbm/science/magnetars/magnl550triggers.html 

13 7 x 1 pixels for XRT WT mode. 
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Figure 1 . Top panel: background subtracted light curve of a burst detected at 02:34:28.194 on 2009 January 22 from SGR J1550— 5418 with the GBM NaI-6 detector. 
It is binned with a time resolution of 8 ms. Middle panel: XRT light curve of the same burst in 8 ms time resolution. Bottom panel: hardness ratios (GBM vs. XRT) of 
the three sub-intervals indicated with the dotted lines. The dot-dashed lines denote the time interval over which the burst spectrum has been accumulated. 


saturated events, in ten other bursts the average readout time 
for central detection areas is smaller than 1 .7 ms, the smallest 
readout time for the WT mode. This indicates that these bursts 
are affected by the photon pile-up in the XRT data. Compared to 
the remaining common events, these bursts have more photons 
in grades higher than 2; as a result an analysis of the data in 
the good grade range (0-2) would lack most information and 
the results would be misleading. Therefore, we also excluded 
these ten piled-up bursts from further investigations. The final 
outcome of all these filters was a selection of 42 bursts, observed 
with both GBM and XRT, that we then used for broadband 
spectral analysis (40 bursts detected on 2009 January 22 and 2 
events on 2009 January 30). 

2.3. Data Reduction 

We extracted time-integrated spectra for our 42 simultane- 
ously detected bursts. In Figure 1, we present the burst detected 
at 02:34:28. 194 on 2009 January 22 as an illustrative example of 
spectral integration ranges. We describe below the procedures 
we followed to generate the GBM and XRT spectra. 

For the GBM data, we selected the Nai detectors with a 
source angle smaller than 60°, and without any blockage from 
the Large Area Telescope or other parts of the satellite, such as 
solar panels and radiators. We determined background levels by 
fitting pre- and post-burst intervals with a first-order polynomial 
using RMFIT \’3.4rcl2 (Mallozzi et al. 2008) and extracted both 
burst and background spectra. We then used grppha to group the 
extracted source spectrum to include at least 15 source counts 
in each energy bin. We generate the response matrices using 
GBMRSP vl.9 for each burst. 

For the XRT data, we selected events with grade range of 0-2 
with xselect, and accumulated the source spectra from a 40 pixel 
long section of the chip centered at the SGR J1550— 5418 lo- 


cation in the same time intervals as used for the GBM spec- 
tra. We extracted background spectra from a region of the 
same size, away from the source. We generated the ancil- 
lary response function file for each burst using xrtmkarf in 
HEASOFT. In our spectral fitting we used the standard re- 
sponse hie “ swxwt0to2s6_20070901v012.rmf ’ provided in the 
Swift calibration database. Finally, we also grouped the source 
spectra to include a minimum of 15 counts in each energy bin. 
We fit all spectra using XSPEC vl2.7.0. 

3. SPECTRAL ANALYSIS 

Motivated by our recent results published in Lin et al. (201 1) 
and van der Horst et al. (2012), we modeled the broadband 
time-integrated spectra of all 42 common events with the 
Comptonized model (COMPT) and the sum of two blackbody 
functions (BB+BB). The COMPT model is in a single power- 
law shape with a high-energy exponential cutoff expressed as 

f = A exp[-£(2 + X)/E pcak ](E/E pn .) 1 , 

where/is the photon flux in photons s -1 cnT 2 kcV ', A is the 
amplitude with units same as f £ pea k is the energy (in keV) at 
which the spectral distribution function peaks, X is the photon 
index, and E p ; v is the pivot energy fixed at 20 keV. The latter 
BB+BB model has been commonly used in the context of SGR 
burst spectra (Feroci et al. 2004; Olive et al. 2004; Israel et al. 
2008). In our spectral fits, we fix the multiplicative interstellar 
absorption term at 3.24 x 10 22 cm -2 , since SGR J1550— 5418 
bursts are short and XRT burst spectra cannot constrain the 
absorption parameter. 14 In our joint fits, we also include a 

14 We adopt here the value obtained with the XRT observations of the 
persistent emission during the burst active episode of SGR J1550— 5418 by 
Scholz & Kaspi (201 1). 
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Figure 2. Spectrum of the SGR J1550— 5418 burst shown in Figure 1. Top two panels show the XRT-GBM joint fit spectrum. Bottom two panels show the GBM-only 
spectrum. The left column are COMPT model fits and the right column are fits with a BB+BB model. The lower parts in each panel show the fit residuals. 

(A color version of this figure is available in the online journal.) 


multiplicative factor to account for the cross-calibration between 
the XRT and GBM spectra. In Figure 2, we show the broadband 
spectral modeling results of the burst of Figure 1 with the 
COMPT and BB+BB models. We note that a BB+BB model fit 
to the combined XRT and GBM data results in better residuals 
than the COMPT ones and are also better than the COMPT fits 
of the GBM data alone (see the two left panels in Figure 2). 
These results indicate that the broadband joint fits can constrain 
better the physical emission model in SGR bursts as we discuss 
in Section 4.1. We calculated the average reduced x 2 for both 
models, and obtained 1.03 for the BB+BB model and 1.11 
for the COMPT model. The standard deviation of the reduced 
X 2 is 0.20 and 0.26 for the BB+BB and the COMPT model, 
respectively. These results further support the BB+BB model 
against COMPT. 

Table 2 presents the joint fit results with Icr errors for all 42 
bursts. Columns 1-3 are the burst numbers, start times in UTC, 
and durations of the time-integrated spectra. The COMPT model 
parameters, i.e., the power-law index (k), E pe and fit statistics 
are shown in Columns 4-6. Columns 7-13 correspond to the 
temperatures of the two blackbody components, the luminosity 
and size of the emitting area for each blackbody component 
(assuming a distance to SGR J1550— 5418 of 5 kpc), and the fit 
statistics of the BB+BB model fits. Since the energy flux values 
obtained with the COMPT model are in agreement with those 


we get from the BB+BB model, we only present in Columns 14 
and 15 the energy flux in the GBM energy band (8-200 keV) 
and the observed flux in the XRT energy band (0.5-10 keV) 
using the BB+BB fits. 

The COMPT model has one less parameter than the BB+BB 
function. Since these two models are not nested, we cannot 
employ a Ax 2 test to determine which model provides a more 
adequate fit. Therefore, we performed extensive simulations 
for each of the 42 bursts. We first fit each joint spectra with 
the seed model (i.e., the model with the smallest resulting x 2 
value): these were the BB+BB model for 33 bursts and the 
COMPT for the remaining nine bursts. Then, we generated 
10,000 simulated spectra based on the seed model of each burst 
and fitted the simulated spectra with both models. For each set 
of 10,000 simulations, we selected fits with well-constrained 
model parameters, requiring the errors of £p ea k of COMPT and 
of the two blackbody temperatures to be less than 20%. We 
constructed the distributions of all fit parameters as well as of 
the fit statistics and fitted them with a Gaussian function. The 
mean values of these distributions for the seed model parameters 
agree with the fit results listed in Table 2. Finally, we calculated 
the percentage of the simulated spectra which result in a smaller 
X 2 value when fitted with the seed model. This percentage, 
defined as the p- value, reflects the significance of the preference 
of the seed model at a given background and fluctuation level. 
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Figure 3. Plot of the p-value vs. the total counts in the GBM energy band 
(8-200 keV). The black dots are the BB+BB bursts, the open circles are the 
intermediate group bursts, and the five-point star is the COMPT burst. 



We conclude that the seed model provides a significantly better 
fit than the other model if p > 0.9. The p - values of all bursts 
are listed in the last column of Table 2. 

We then grouped the 42 bursts into three categories 
based on their resulting /7-values: the BB+BB burst, where 
BB+BB is significantly preferred to the COMPT model; the 
COMPT burst, where the COMPT model is significantly 
better than BB+BB; and the intermediate group containing 
bursts for which both the BB+BB and COMPT models pro- 
vide equally acceptable fit results ( p < 0.9). We find that 
31 events are the BB+BB bursts, only one is a COMPT 
burst (event 18), and 10 are in the intermediate group. Fig- 
ure 3 displays the /7-values of all bursts versus their total 
GBM counts. Note that a similar trend is obtained if the 
p - values were plotted with their corresponding total XRT counts. 
We find that the bright bursts prefer the BB+BB model to the 
COMPT model. This might be an indication of higher opacity, 
on average, in the more luminous bursts. We explore both the 
statistical and the correlative behavior of all model parameters 
in the next sections. 

Our joint spectral fits provided the opportunity to investigate 
the cross-calibration of the XRT and GBM instruments with 


the two spectral models. We determined a multiplicative factor 
between the XRT and GBM detectors for each model. The left 
panel of Figure 4 shows the behavior of this factor from the 
BB+BB model fits as a function of total counts in the brightest 
GBM detector. The same plot for the COMPT model fits is 
presented in the right panel of Figure 4. The instruments are 
perfectly cross-calibrated if the multiplicative factor is equal to 
1 (the dotted lines in Figure 4). For brighter bursts, the factor 
from the BB+BB model fit is better constrained, and it does not 
change significantly from burst to burst: its weighted mean value 
and I a error are 1 . 17 ±0.05. However, the weighted mean of the 
COMPT model derived factor is 0.58 ± 0.03, much smaller and 
further from 1 . We note that the values obtained from a better fit 
model are closer to the perfect cross-calibration factor of 1 . We 
find that the constant factor is not significantly correlated with 
any spectral parameters. We conclude that the cross-calibration 
of the XRT and GBM instruments works reasonably well within 
the fluence and energy range of the bursts in our sample. Future 
efforts to better understand the cross-calibration of these two 
instruments should include bursts with wider energy and fluence 
ranges, and perhaps, different spectral models. 

3.1. Comptonized Model 

We present the distribution of power-law indices obtained 
from the joint fit of the XRT and GBM spectra (histograms with 
thick lines) in the left panel of Figure 5. We also present in the 
same figure the distribution of the indices as obtained by fitting 
the GBM spectra only (histograms with thin lines). It is clear 
that the latter fits yield, on average, lower power-law indices 
(i.e., the spectrum is harder), reflecting the overall broadband 
curvature of the burst spectra. We fit each distribution with a 
normal function and find that the joint broadband fit mean index 
value is —0.58 ± 0.09 (width of 0.43 ± 0.11, dotted lines in 
Figure 5, left panel), while the GBM only fit mean index is 
—0.87 ± 0.05 (width of 0.42 ± 0.06, dashed lines in Figure 5, 
left panel). We show the distribution of the £p e ak values obtained 
from the joint fit in the right panel of Figure 5 . A normal function 
fit to this distribution yields a mean of 45.0 ±2.1 keV with a 
width of 10.9 ± 2.2 keV. The £p e ak values from joint fit agree 
with those from GBM data only fit very well, see detailed 
discussion in Section 4. 1 . 

In Figure 6, we present a plot of the joint fits £p e ak values 
versus the observed fluence/flux of bursts in the 0.5-200 keV 




Figure 4. Plot of the multiplicative factor values from the BB+BB model (left) and the COMPT model (right) fits vs. the total counts (8-200 keV) in the brightest 
GBM detector. The black dots are the BB+BB bursts, the open circles are the intermediate group bursts, and the five-point star is the COMPT burst. The dotted line 
indicates the factor being equal to 1, and the solid and dashed lines display the weighted average of the factor and its la error, respectively. 


7 


The Astrophysical Journal, 756:54 (12pp), 2012 September 1 


Lin et al. 



Index X E peok (keV) 


Figure 5. Distributions of the COMPT model index (left) and Fpeak (right). The dotted lines show the best-fit normal functions. The thick histograms in the left panel 
are the index distribution of the XRT-GBM joint fit, while the thinner histograms are the same distribution for the GBM data only fit. 



io ' 7 icr 6 



Energy Fluence (erg/cm 2 ) Energy Flux (erg/cm 2 /s) 

Figure 6. Plot of Fpeak vs. observed energy fluence (left) and average flux (right) between 0.5 and 200 keV. The black dots are the BB+BB bursts, the open circles are 
the intermediate group bursts, and the five-point star is the COMPT burst. 


range. We do not find any anti-correlation or a broken power-law 
trend as seen by fitting the GBM data of SGR J1550— 5418 (van 
der Horst et al. 2012) and SGR 0501+4516 (Lin et al. 2011), 
which we attribute to the fact that the 42 common events cover a 
much narrower fluence range, about 1/3 of that in the complete 
GBM burst sample from the same active period (van der Horst 
et al. 2012). 

3.2. Two Blackbody Model 

From simulation, we find that the BB+BB model provides 
a significantly better fit to 31 common events. We use these 
BB+BB bursts to investigate the properties of the model param- 
eters. In Figure 7, we present the distribution of temperatures of 
the cool (left panel) and hot (right panel) BB components, re- 
spectively. We fit a normal function to these distributions and ob- 
tain a mean value for the cool BB temperatures of 4.4 ± 0.2 keV 
(width, 0.8 ±0.1 keV), and for the hot BB temperatures of 
16.0 ± 0.4 keV (width 2.2 ± 0.4 keV). We also calculated the 
weighted mean values for cool and hot BB temperatures, which 
are 4.2 ±0.1 keV and 14.8 ± 0.2 keV, as well as their standard 
deviations, 0.9 keV and 2.7 keV. These mean temperatures are 
in agreement within uncertainties with the values obtained by 


fitting only GBM spectra of SGR J1550— 5418 bursts (van der 
Horst et al. 2012) and similar to the BB+BB temperature values 
obtained for the bursts of other SGR sources (Feroci et al. 2004; 
Olive et al. 2004; Israel et al. 2008; Lin et al. 201 1). 

In Figure 8, we demonstrate the anti-correlation between 
the emission area and the temperature of the cool and hot 
BB components. The Spearman rank order correlation test 
yields a correlation coefficient for the hot BB of —0.79 with 
a chance probability of 1.40 x 1 0 7 . The correlation for 
the cool BB is not as significant as the hot component: its 
Spearman correlation coefficient is —0.62 with a probability 
of 1.75 x 10~ 4 . We fit the emission area versus temperature 
for the cool and hot BB with power laws, and obtain power- 
law indices of —1.5 ± 1.4 and —4.5 ± 0.9, for the cool and 
hot components, respectively. We also fit the emitting area 
versus the two BB temperatures together with a single power 
law and obtain the best-fit power-law index of —3.5 ± 0.2 
(shown as the solid line in Figure 8). Note that this value 
is very close to the theoretical expectation from a single BB 
with fixed luminosity, R 2 oc ( kT)~ 4 . The observed departure at 
the high-temperature end from this ideal form reflects the higher 
luminosity present in the hot BB component relative to the cooler 
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Figure 7. Distributions of the temperatures of the two blackbody components of the BB+BB model. The dotted lines show the best-fit normal functions. 
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Figure 8. Emission area as a function of blackbody temperature for both 
blackbody components of the BB+BB model fit. The hot and cool blackbody 
components are displayed with triangles and circles, respectively. The filled 
symbols are the BB+BB bursts. The intermediate group bursts are shown as 
open symbols. The solid line indicates the R 2 oc (kT) 3,5 relation, the best-fit 
power-law function with all emission areas and temperatures for BB+BB group 
bursts. 
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Figure 9. Correlation between the total energy emitted from the hot and cool 
blackbody components. The filled and open circles present the BB+BB and 
intermediate group bursts, respectively. The best power-law fit (E hot oc T.ff ) 
for the BB+BB group is shown as a solid line. The dotted line is where the two 
BB components have equal energy. 


one; on average, the hot BB energy is about twice the one emitted 
from the cool BB. The two energies are highly correlated with 
a Spearman correlation coefficient of 0.88, corresponding to a 
chance probability of 6.45 x 1 0“ 1 1 ; a power-law fit yields an 
index of 0.99 ± 0.05 (Figure 9). 

4. DISCUSSION 
4.1. Impact of the XRT Data 

We have used two different instruments, the Swift/X RT 
and the Fermi /GBM, to perform a time-integrated broadband 
spectral analysis of 42 common events from SGR J1550— 5418. 
By adding the XRT data, we extended the lower energy bound 
of our earlier spectral analysis (van der Horst et al. 2012) from 
~10 keV to ~0.5 keV. For most model parameters, our joint fit 
results agree well with the results from fitting the GBM spectra 
only, as shown in Figure 10. The correlation coefficients between 
all but one (the COMPT power-law index) model parameters 
derived with and without the XRT data are larger than 0.94, 
corresponding to a probability smaller than 3.3 x 10” 1() . The 


average COMPT index without the XRT data is —0.87 ± 0.05, 
consistent within errors with the mean of 0.92 obtained 
from a much larger burst sample (van der Horst et al. 2012). 
However, the inclusion of the XRT data better constrains the 
COMPT indices, which become harder than the ones derived 
from spectral fits to the GBM data alone, as shown in Figure 5. 
Therefore, we conclude that the COMPT fit to the GBM data 
only overestimates the emission in the lower energy bands. 

Our analysis provides an important diagnostic for the model 
preference between the COMPT and BB+BB models. By adding 
the XRT data, we find that 31 of 42 bursts are statistically better 
described by BB+BB. This fact, combined with the observed 
steepening of COMPT model indices when excluding XRT data, 
highlight the generic broad curvature of the SGR J1550— 5418 
burst spectra. We note here that the joint analysis of XRT 
and GBM spectra is limited to an absorbed energy fluence 
range of 8.3 x 1CT 8 -1.5 x 10“ 6 erg cm" 2 (0.5-200 keV), 
since the brighter events would cause pile-up in the XRT data 
and the dimmer events would not yield high enough statistics 
in the XRT data for a constraining spectral analysis. We discuss 
the theoretical implications of these results in the next session. 
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Figure 10. Plots of correlations between model parameters obtained by the joint XRT-GBM fits and the parameters obtained by fitting the GBM data only. The black 
dots are the BB+BB bursts, the open circles are the intermediate group bursts, and the five-point star is the COMPT burst. The dotted lines represent the x = y trend. 
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Figure 11. Top: phase probability density profile of the BB+BB group bursts from SGR J1550— 5418 ; middle: count phase distribution; bottom: persistent emission 
pulse profile from contemporaneous XMM observations in the 0.5-10 keV band. The horizontal dashed lines in each panel represent the mean value of the burst phase 
probability density (top panel), the burst phase counts (middle panel), and the persistent emission count rate (bottom panel). 


We also noted the fact that Israel et al. (2008) investigated 
broadband spectral properties of a very rare event (the storm) 
from SGR 1900+14, while our investigations are about much 
more common typical short bursts. Therefore, our results extend 
the Israel et al. (2008) results to the more common magnetar 
outbursts. 

We investigated here SGR J1550— 5418 burst spectra in a 
time-integrated manner. As seen in the bottom panel of Figure 1 , 
the hardness ratios of the three parts of the burst, designated by 
the three peaks in the XRT and GBM light curves, show a clear 
hard to soft spectral evolution. Detailed time-resolved spectral 
analysis of SGR bursts would provide important insight to the 
spectral evolution of SGR bursts with time. 


4.2. The BB+BB Model 

The most plausible interpretation of the BB+BB model is the 
emission originating from two hot spots with different tempera- 
tures near or on the neutron star surface or in its magnetosphere 
where local thermodynamic equilibria are achieved. It must be 
emphasized that a BB spectral fit is an idealization to that emit- 
ted by the physical environment of a real photosphere. Due 
to possible gradients of temperature with optical depth into an 
evolving region that is approximately in local thermodynamic 
equilibrium, significant distortions from a true blackbody form 
are predicted in SGR photospheric spectral models (e.g., Ulmer 
1994; Thompson & Duncan 1995; Lyubarsky 2002; Israel et al. 
2008). 


To better understand the BB+BB behavior and uncover its 
relation with the spin properties of SGR J1550— 5418, we 
investigated the phase characteristics of the 3 1 BB+BB bursts, 
as follows: We first selected all XRT counts collected during 
31 burst intervals and converted their arrival times from the 
Swift mission time to the corresponding time at the solar system 
barycenter. We then calculated the spin phase for each burst 
count using the appropriate spin ephemeris of epoch (MJD) 
54854 as reported by Dib et al. (2012) using both RXTE and Swift 
observations. We present the phase distribution of burst counts 
in the middle panel of Figure 1 1 . To ensure that the distribution 
is not dominated by the excessive counts of the brightest bursts, 
we also calculated the probability density for each phase bin, 
which is the average of the normalized (by total counts) phase 
distributions for all bursts, as shown in the top panel of Figure 1 1 . 
We find that the probability distribution of the burst counts is 
not uniform over the spin phase of SGR J1550— 5418 and the 
deviation from the mean probability is significant: we calculate 
the root-mean-square deviation of the phase probability density 
function from its mean as 0.021 ± 0.001. We also compared the 
phase probability density function to the persistent emission 
phase profile (bottom panel in Figure 11) obtained using 
contemporaneous XMM observations (Dib et al. 2012). The 
phase probability density function is marginally anti-correlated 
with the persistent emission phase profile in our burst sample 
with the correlation factor of —0.5 corresponding to a chance 
probability of 3.4 x 1 0 : . This indicates that the burst emission 
regions on the neutron star surface are not necessarily associated 
with the site persistently emitting in X-rays (typically a BB with 
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a temperature of 0.5 keV). This is in agreement with the crustal 
fracturing mechanism for SGR bursts (Thompson & Duncan 
1995; Braithwaite & Spruit 2006; Perna & Pons 2011) as any 
portion of the solid crust can fracture if the magnetic stress built 
up is near the threshold to rupture. We also find that the burst 
probability of some spin phases in SGR J1550— 5418 is higher. 
This could be attributed to a non-uniform surface magnetic field, 
with some regions having larger magnetic stresses than others. 

4.3. The COMPT Model 

A Comptonization spectrum emerges when low-energy pho- 
tons are repeatedly upscattered by the thermal electrons in a 
corona until the photon energy reaches E ~ kT e . We find that 
the mean value of the Comptonization peak energy, £p e ak, is 
44.8 keV, which indicates an average temperature of the ther- 
mal electrons, (T e ) ~ 5.1 x 10 s K. In other words, the average 
speed of the electrons in the corona is ~0.4 c, where c is the 
speed of light. 

The essence of Comptonization spectra is discussed at some 
length in Lin et al. (201 1), in the context of GBM observations 
of SGR J0501+4516 bursts. While the turnover energy pro- 
vides a diagnostic on the hot electron temperature, the power- 
law slope below the vF v peak energy defines a measure of 
the opacity in the Comptonizing region. Specifically, in the sim- 
plest theoretical constructs (e.g., see Rybicki & Lightman 1979) 
the index X for a differential photon spectrum dN/dE oc E x 
depends only on the magnetic Compton y-parameter y B = 
4kT e /(m e c 2 ) max{r B , r^} via X = 1/2 — ^9/4 + 4 /y B . Here 
max{r, r 2 } is the mean number of scatterings per photon by the 
hot electrons, where r is the effective optical depth for scat- 
tering, which in our case is modified by the strong magnetic 
field and thus dubbed the magnetic optical depth and denoted 
by x B . Since we require the effective value of r«, one needs 
to calculate some Rosseland-type mean opacity, averaged over 
photon angles and polarizations and accounting for the effects 
of the strong magnetic fields present, such as anisotropy and 
polarization-mode switching through scattering in non-uniform 
B (see also the discussion in Lin et al. 2011). This index is 
realized only in the energy range somewhat above the soft 
photon injection energy, E x , presumed to be surface thermal 
X-rays, and somewhat below the characteristic energy associ- 
ated with the hot thermal electrons, in this case marked via the 
vF v peak energy £p ea k- Moreover, the above relationship for y B 
requires the scattering to be in the Thomson regime, and the 
mean photon energy to be lower than that of the electrons. 

It is evident that y B 1 cases yield the flattest Comptonized 

spectra with index around X 1. The COMPT fit indices in 

Table 2 as well as its distribution in the left panel of Figure 5 
are nearly always harder than this, indicating that repeated 
Compton upscattering has difficulty in generating the observed 
flat spectra. We note parenthetically that this was also the case 
for around 1/3-1/2 the bursts reported for SGR J0501+4516 
in Lin et al. (201 1). The inclusion of XRT data in the current 
work extends the spectral coverage to much lower energies thus 
enabling us to determine the value of X significantly better 
than in previous works. This results in systematically higher 
values, corresponding to a harder spectral slope, with a mean 
value of —0.6 ±0.10 instead of —0.92 ± 0.05, which has made 
this problem worse. The fact that a large y B is demanded in 
this fitting protocol would suggest high opacity and strong 
thermalization might be active in the burst emission region. 
In such cases, the above dependence of y B on T e and r B is not 
operable, and the Comptonization is saturated, and described 


instead by a modified Wien or modified BB spectrum (e.g., see 
Rybicki & Lightman 1979). Moreover, the upscattered power- 
law photon population is generally substantially lower in total 
number relative to the seed surface X-ray population (i.e., it 
presents itself as an X-ray tail), a broadband spectral shape 
that is at odds with the spectral curvature inferred from the fits 
here. Hence, there is no strong mandate to prefer a classical, 
unsaturated Comptonization model for the bursts reported 
here. Accordingly, the spectra might naturally be expected to 
exhibit more truly thermalized character, for example, a two- 
blackbody or multi-blackbody signal emanating from ar B » I 
zone, possibly eliciting spectral distortion imposed by transport 
within the photosphere (e.g., see Ulmer 1994; Lyubarsky 2002). 
Perhaps this is what this broadband XRT/GBM analysis has 
enabled for a magnetar for the first time: the clear discrimination 
between COMPT and BB+BB spectral models. 
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